Hidden Symmetries of Electronic Transport in a Disordered One-Dimensional Lattice 
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Correlated, or extended, impurities play an important role in the transport properties of dirty 
metals. Here, we examine, in the framework of a tight-binding lattice, the transmission of a single 
electron through an array of correlated impurities. In particular we show that particles transmit 
through an impurity array in identical fashion, regardless of the direction of transversal. The 
demonstration of this fact is straightforward in the continuum limit, but requires a detailed proof 
for the discrete lattice. We also briefly demonstrate and discuss the time evolution of these scattering 
states, to delineate regions (in time and space) where the aforementioned symmetry is violated. 

PACS numbers: 72.10.Bg, 72.10.Fk 
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I. INTRODUCTION 

In a one-dimensional lattice, correlated impurities or 
defects may give rise to extended electronic states for par- 
ticular energy levels Various disordered systems with 
correlation have been investigated, such as the random 
dimer^i the random trimer^ non-symmetric dimersp 
the random trimer-dimerj& the Thue-Morse lattice^ and 
the random polymer chain£ A built-in internal struc- 
ture of the impurity configuration is important*^ even if 
an internal symmetry is not necessary^ In the context of 
electronic transport, the transmission resonance for a cer- 
tain energy value implies the disordered system behaves 
like an ordered lattice for the corresponding electronic 
state, which is extended throughout the entire lattice. In 
this respect, it is possible to understand the increase in 
the conductivity for a conducting polymer 2 . 

One may, in fact, envision many other impurity con- 
figurations that give rise to extended states for various 
energy values. For example, let us consider a pair of im- 
purities with a potential V with respect to the lattice 
potential (taken to be zero). We assume that the spac- 
ing between the two impurities is always n^a, where a 
is the lattice constant and rid is a positive integer. This 
is then a 'dimer', whose length can be any value. If we 
adopt a tight-binding model, so that energies are given 
by E — — 2iocos(fca), where to is the nearest-neighbor 
hopping amplitude, then one can obtain the transmis- 
sion probability \T\ 2 for this state. It is 
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cos(fcn d ) + f^ 
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where to and a have been set to unity for simplicity. 
Inspection of Eq. JIJ reveals a critical strength of the 
impurity potential V c /to = ±2 /rid for repulsive and at- 
tractive interactions, respectively, which determines how 
many states are extended, i.e. have unit transmission. A 
graphical construction readily shows that, for V > 2/n<j, 
there are rid — 1 extended states while for V < 2 /rid, there 
are rid extended states. As in the nearest-neighbor dimer 




FIG. 1: (Color online) Transmission probability as a function 
of E for randomly distributed 50 impurity dimers (blue) and 
500 dimers (green) in comparison with a single dimer (red). 
The inner spacing of a pair na — 3 and the potential V = 1/2 
while the critical value Vo = 2/3. The transmission proba- 
bility is calculated as an ensemble average over 50 impurity 
configurations. 



case, consideration of one 'dimer' helps us to understand 
electronic transport in a lattice with many 'dimers', fol- 
lowing the analysis of Ref. 1 . In Fig. 1, we show |T| 2 
as a function of E for randomly distributed 50 impurity 
dimers (blue) and 500 dimers (green) in comparison with 
a single dimer (red) . The inner spacing of a pair rid — 3 
and the potential V = 1/2. Note that in this instance 
V < V c = 2/3. The transmission probability is calcu- 
lated to be an ensemble average over 50 impurity con- 
figurations. Remarkably, the transmission remains unity 
at the three energies for which a single 'dimer' has unit 
transmission, and the 'width' of the resonance narrows 
as the number of impurities increases. 

This specific example is given by way of introduction to 
more general considerations of impurity potentials on a 
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lattice. In this paper we will show there exist hidden sym- 
metries underlying the physics of electronic transport in 
a one-dimensional lattice with an arbitrary impurity con- 
figuration, using the transfer matrix formalism. We will 
also demonstrate the time evolution of a wave packet by 
direct diagonalization of the Hamiltonian since the trans- 
fer matrix formalism does not produce the dynamics of a 
wave function. While the dynamics can be insightful, we 
are also hopeful that present-day technology can spatially 
and temporally resolve some of the interesting dynamics 
that result from our (and future) work. 

Perhaps the most general correlated impurity config- 
uration is the following. Suppose there are N impuri- 
ties AiA 2 • • -An embedded at the sites 1 through N in a 
host lattice. The corresponding impurity potentials are 
given by (U\, U 2) • • • , U N ), where Ui (i = 1, • • • , N) 
can be positive, zero, or negative. Our main result is 
that the transmission amplitude T for this impurity con- 
figuration is identical to T for its reverse configuration; 
namely, T(U U U 2 , U N )=T(U N , U N _ U Pi). 
Given this then obviously the equality holds for the trans- 
mission |T| 2 and reflection |P| 2 probabilities. In other 
words, |T| 2 and |P| 2 are independent of the direction of 
the incoming wave. This symmetry is striking because 
the impurity configuration is assumed arbitrary without 
possessing any parity. For example, it has been shown 
that the chains of ABBABAAB- ■ • and BAABABBA- ■ • 
contribute identically for electronic transport 7 . This 
is a simple example of the general symmetry we de- 
scribe. Note that an exchange of any two impurity 
potentials does not show this symmetry; for example, 
T(U U U 2 , U 3 , U 4 ) ? T(U 2 , U u U 3 , P 4 )2. 

We should emphasize now that this symmetry does 
not hold in the impurity region (as perhaps one would 
expect), but this can only be determined by solv- 
ing the time-dependent Schrodinger equation (see be- 
low). We also explain another symmetry associated 
with a particular momentum k = n/(2a), where a 
is the lattice constant: \T(Ui, U 2 , PaOI 2 = 

\T(—Ui, —U 2 , — Un)\ 2 - This equality indicates a 

potential barrier or well produces identical transmission 
and reflection probabilities for the electronic state with 
k = n /(2a). This is unrelated to the other remarkable 
feature of wave packets on a one-dimensional lattice with 
nearest- neighbor hopping and k = n/(2a): they do not 
diffuse as a function of time, because they lack group ve- 
locity dispersion. Indeed, one can show that no matter 
what the range of the hopping is, there will be a wave 
vector for which the wave packet retains its initial width. 



II. FORMALISM 

We start with a tight-binding Hamiltonian for a one- 
dimensional lattice with N impurities: 



where the hopping amplitude to will be set to be unity, 
cf creates an electron at a site i, and X represents an 
impurity configuration. The Schrodinger equation is 
H\ip) — E\ip) with \tp) = J2 3 V'jCjNO}; this becomes 



[ip j+1 + tpj-x] + Ujipj = Eijjj 



(3) 



where tpj is an amplitude to find an electron at site j. 
In order to use the transfer matrix formalism we write 
Eq. Q in a matrix form as follows: 



ipj+i 

i>3 



U t -E -1 
1 



tpj-1 



Mi 



(4) 

Note that Mj is a unimodular matrix, i.e. det(Mj) = 1. 
The wave functions t/jl (for i < 1) and ipR (for i > N) are 
ipL = e lkXi + R e~ lkXi and ip R = T e lkXi , where Xi — a-i. 
Using the transfer matrix formalism, one can express the 
coefficients R and T in terms of k, Ui, and E as follows: 



T 
iT 



= P 



1 + R 

i(l - R) 



where P = S^MS with S = 



cos(fc) sin(fe) 



(5) 



and 



1 

M = MatMjv-1 • • • Mi. It is misleading to express M = 
HfL 1 Mi because Mj and Mj (i ^ j) do not commute 
with each other. Solving Eq. JSJ, one can obtainiS 



T 
R 



2i 



l(Pll+P 2 2)+Pl2-P21 

Pi2 + P21-1 (Pn - P22) 

% (Pn + P 22 ) + P 12 - P 2 i 



(6) 
(7) 



Consequently, the transmission probability |T| 



4/ 



tr(PP) 



and the reflection probability \R\ 



H 



-tt 



tr(PP) - 2 / tr(PP) + 2 , where tr(P) means the 

trace of P, and P is the transpose of P. Note that 
Eq. is readily obtained with this formalism by us- 
ing M' = MvMq^Mv, where My = M X (U X = V) and 
Mq = Mi (Pi = 0), and is the inner spacing of the 
impurity pair. 

One of the symmetries we introduced earlier is 
T(Ui, U 2 ,---, U N ) - T(U N , Pat-!,---, Ui). In or- 
der to show this equality, let us introduce Q = S~ 1 M'S 
where M 1 = M X M 2 ■ ■ ■ M N . Notice the order of the ma- 
trix multiplication because Mi and Mj (i ^ j) do not 
commute with each other. Since the desired equality im- 
plies that T[P] = T[Q], we need to show Pn + P 22 = 
Q11+Q22 andP i2 -P2i = Qi 2 — Q 2 i- On the other hand, 
from the definitions of P = S^MS and Q = S^M'S, 
this equality in turn implies that Mn = M n , M 22 — 
M22, and M12 - M 2i = M' X2 - M 21 . Note, however, that 
in general M ^ M', as we will show later. 

For simplicity, we define Ui — Ui — E. 



[c+Cj+i +cf +1 Ci] +22u i cfc i (2) troduce 2x2 matrices a and such as a = 



We also in- 
'1 0\ 
j 



3 



and j3 = 



-1 

1 



to have M, = UiCt + (3. Note that 



M = A/jvMtv-i ■ ■ • Mi, we obtain 



a 



a, /? 2 



-1, a(3 + /3a = /?, and a(3a = 0. Since 



M = {U N a + I3){U N - I a + (3)---(U ia + !3) 

N N 

= P N + YY ( U n u n ■ ■ -UjJ^-^a^-^-^a ■ ■ ■ afi^-^ufi^ 1 
n {j„} 



(8) 



where {j n } means ji,j 2 , ■ ■ ■ , j n = 1)2, , • • • , N with ji < j 2 < ■ ■ ■ < j n - Since a(3a — and [3 2 = — 1, non- vanishing 
terms should have ji — — 1 = even, where I — 2, 3, • • • n for a given n. Therefore, we obtain 



N N 



M = /3 N + Y Y (UhUj2---Ujn)(- 1 )^ jn ~ h ~ n+1) P N ~ jn aP jl ~ 1 



0) 



r 



where {{jn}} means {j n } with ji —ji-i — 1 = even. Since 
1 0\ „ f -l\ n f Q 





-1 



1 



and 



, we use these as the basis matrices to 



expand M: 



M = c\a + c 2 a(3 + cspa + c±(ia.f3 



c\ -c 2 

C 3 -C4 



(10) 



As we mentioned, the equality of the transmission am- 
plitude, T[P] = T[Q], is associated with some particu- 
lar relations among the components of the two matrices 
M = M\M 2 ■ ■ ■ M N and M' = M N M N - X ■ ■ ■ M x . Follow- 
ing the same way as for M, we expand M 1 in terms of 
Ui (i = 1, 2, • • • , N), a, and {3 to obtain 



M' = {U ia + f3){U 2 a + f3)---{U N a + 13) 

N N 

= P N + YY ( u h u h ■ • •«i„)/P 1 ~ 1 a0 fe ~ il-1 a ■ ■ ■ a p in - Sn - 1 - l ctP N - jn . 
It is also true for M 1 that non-vanishing terms have ji — ji-i — 1 = even because a/3a = 0. Therefore, 

N N 
" {{jn}} 



(11) 



(12) 



r 



Expanding M', again, in terms of the basis matrices, we To obtain i) c\ of M, we should have jt — 1 = even and 

N — j n — even. Similarly, ii) for c 2 , ji — 1 = odd and 
N — j n — even, iii) for C3, j\ — 1 = even and N — j n = 
odd, and iv) for C4, — 1 = odd and N — j n = odd. In 



know 

M' = cia + c 2 a/3 + c^a + c'^aP 



c l c 2 



(13) 



Comparing Eqs. (|10Jl and 113|l . the equalities Mn = M( 1 , 
A/22 = M22, and M12 — M21 = M[ 2 — M 21 are equivalent 
to the equalities c\ = c' l7 C4 = c' Al and C2 + C3 = c 2 + c 3 . 
Since (3 N does not depend on Ui and /3 W = (-l) w / 2 or 
(— l)^^ 1 )/ 2 /? for even (or odd) TV, we can ignore f3 N for 
the purpose of showing these equalities, or, alternatively, 
we can redefine M as M - f3 N and M' and M' — (3 N . 



terms ofUi, we obtain 

N N 

n {{i n }}i 

iV JV 

C2 = Y, Y ( U h U h 

« {{Jn}}2 



■Z4J(-i)' (JV_n) 
•^J(-i)^ (Ar ^ 1) 



4 



C"3 



(JV-ra-1) 



TV TV 

E 

" {{>Jb 

TV TV 

C4 = E E (^ 2 --^j(-i)i^-™- 2 ) 

where {{j„}}i means {{j„}} with {j x - 1,N - j n ) = 
(even, even). Similarly, {{j, l }} 2 = {{jn}} with (odd, 
even), {{j„}} 3 = {{jn}} with (even, odd), and {{j„}} 4 = 
{{in}} with (odd, odd). 

For c'i, c' 2 , c 3 , and c 4 of Eq. I|13[) . we should have 
(ji — l,N — j n ) — (even, even), (even, odd), (odd, even), 
and (odd, odd), respectively. Note that we need (even, 
odd) for c' 2 while (odd, even) for c 2 . On the other hand, 
we need (odd, even) for c 3 while (even, odd) for C3. 
Therefore, we obtain 



TV 



TV 



I 

c l 



£ £ (^ v ..^)(-i)^«) 

« {{jn}h 



(u jl u h --.u jn )(-i)^ N - n -^ 



I 

c 4 



TV TV 

E E 

" {{jn}h 

TV TV 

£ £ (^-•^•J(-1)' (JV — X) 

™ {{J,.}}2 

TV TV 

E E (%%-z4J(-i) i(JV — 2) 

™ {{j„}}4 



Consequently, we have c\ — c[, C4 = C4, and C2 + C3 = 
c' 2 + c 3 ; in other words, T[P] = T[Q], or equivalently, 
T(U U U 2 , C/tv) = T(U N , C/tv-i, E/i) In fact, 
M12 = — M 21 and M 2 \ = —M[ 2 . Note, however, that in 
general M ^ M' . 

One may anticipate a similar symmetry in the con- 
tinuum limit iii Consider a time-dependent Schodinger 
equation: idtip(t,x) — H(x)if)(t,x), where the Hamil- 
tonian H(x) includes an impurity potential V(x) for 
x < I. The wave function for |x| > I can be rep- 
resented as V(M > = [e lkx + Re- lkx ] Q(-x - I) + 
Te lkx Q(x — I), where 0(x) is the step function. Let 
us introduce another time-dependent Schodinger equa- 
tion: id t ip'(t,x) = H(—x)ip'(t,x), where ip(\x\ > I) — 
[e lkx + R' e~ lkx ] B(-x - I) +T' e tkx @(x - I). Invok- 
ing both space inversion and time reversal, one can see 
that ip(t, —x) and ip*(—t, —x) satisfy the same equation 
as ip'(t,x) does. Equating aiip(t, —x) + a 2 ip*(—t, —x) — 
f(t, x) for \x\ > I leads to T 1 = T and R'T* + R*T = 0. 
The first relation corresponds to T[P] = T[Q] in a lattice. 
The second relation, however, does not hold in a lattice. 
Instead, we found in a lattice 



T*[P] (R[Q]-R[P}) 



n 2ik 



T[P] {R* 



R*[P}) 



(14) 



Nevertheless, we can still define R[Q] = e iS R[P}. The 
phase S is determined by e lS — -e 2lk R*T/(RT*) in a 



We also found another symmetry associated with elec- 
tronic transport in a one-dimensional lattice; in this case 
there is no applicability in the continuum limit. Suppose 

k = tt/2; then 5 = (® and 5 = 5 = S~K Note 

that for a symmetric (2 x 2) matrix X such that X = X , 
Mi satisfies MiXMi = MjXMj with Ui -> -Hi. Now let 
us consider tr(PP): 

tr (pP) = tr (m n ■ ■ ■ Mi Mi • • • Mjv 
= tr ( M N ■ ■ ■ M X M X ■ ■ ■ M N 



tr (QQ) = tr {QQ)_ u 



(15) 



This symmetry indicates that for k — tt/2, \T\ 2 and |i?| 2 
depend only on \Ui\. If E = — 2cos(fc), then Ut = Ui. In 
this instance, |T(Z7i, • • • , U N )\ 2 = \T(-U N , ■ ■ ■ , -U x )\ 2 . 
Since |T(Ui, ■ • ■ , U N )\ 2 = \T(U N , ■■■,U 1 )\ 2 from a gen- 
eral derivation, \T(Ui, • • • , U N )\ 2 = \T(-U h • • • , ~U N )\ 2 . 
Consequently, a potential barrier and a potential well 
give identical transmission and reflection probabilities for 
k = tt/2 and E = —2 cos(fc). 



III. TIME EVOLUTION OF A WAVE PACKET 

So far we have been using the transfer matrix for- 
malism. Since, however, the transfer matrix formalism 
is based on the time-independent Schrodinger equation, 
we cannot see the time evolution of a wave function, 
which may indicate the significance of the symmetries 
we have shown. Let us consider a wave packet with aver- 
age position xq and average momentum ko at time t = 0: 
l*(0)) = Ei^i,0)c+|0), where 



(p(xi,0) 



(27ra 2 )V4 



^ik (xi-xo) -j(xi-xo) 2 /a 2 



(16) 



We introduced the initial uncertainty a associated with 
position. We wish to propagate (in real time) this wave 
packet towards the potentials. To do this we first diago- 
nalizc the Hamiltonian, which is an (L x L) matrix, where 
L is the total number of lattice sites under consideration. 
We thus obtain the eigenstates |n) and the corresponding 
eigenvalues e„ such that H\n) — e n \n). An eigenstate is 
a column vector with L components which describes the 
probability of finding an electron at a particular site in 
the lattice. With eigenstates in hand, the time evolution 
of the wave packet is given by 



|*(t))=£>><n|*(0)>e- 



ie n t 



(17) 



lattice while e 



,6 



-R*T / (RT*) in the continuum limit. 



As time goes on, the wave packet initially at xq moves 
to the potential region and scatters off the potential. A 
part of the wave packet is reflected and the other part 
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Fig.2 ( Kim et al ) 




depending on the impurity configurations. In particular, 
the scattering completely differentiates [I] from [II], as 
indicated by the very differently behaved red and blue 
curves in the scattering region. Nonetheless the proba- 
bility that emerges (either in transmission or reflection) 
is identical for both, in agreement with the symmetry 
we just proved. Note that the transmitted wave packets 
are identical in all other respects as well whereas the re- 
flected wave packets have identical shape, but are phase- 
shifted with respect to one another. In fact, one can 
show that the phase shift originates from the phase 5 in 
R[Q] = e lS R[P\. The shift measured by the difference 
between the two reflected wave packets for [I] and [II] at 
their half width is determined by dS(k)/dk\ k . 

IV. CONCLUSIONS 



FIG. 2: Time evolution of a wave packet in two impurity 
configurations [I] (red solid) and [II] (blue dashed curve). 
Each configuration has 5 impurities: for [I] the potential is 
(—0.5, 0, 0.5, 1.5, 1) while for [II] its reverse holds. 



is transmitted. Mathematically we can define the re- 
flection and the transmission as \R\ 2 = J2i e ft \<f{ x iit)\ 2 
and |T| 2 = Y, r i g ht\ L P^ x ^ t )\ 2 ^ respectively as t -> oo, 
where Y^uftCright) means the summation includes only 
the left(right) side of the impurity region. 

Fig. 2 shows the time evolution of a wave packet 
impinging on 5 impurities located at sites 300 to 304. 
The parameters of the wave packet are given as follows: 
Xq = 150, a = 20, and ko = 7r/3. Note we choose a — 20 
so that the size of the wave packet is much greater than 
the size of the impurity region. This means that we can 
consider the wave packet as a plane wave in this case. 
In fact, the numerical results of the transmission and re- 
flection probabilities are in close agreement with those 
obtained through the transfer matrix analysis. The im- 
purity configuration is determine by 5 impurity poten- 
tials: (Ui, U 2 , U 3 , J7 4 , U 5 ). For [I] (red solid), we 
have (-0.5, 0, 0.5, 1.5, 1) while for [II] (blue dashed 
curve), we have the reverse order. As clearly shown in 
Fig. 2, the time evolution of the wave packet is different 



We have determined the transmission and reflection 
characteristics for various impurity configurations on a 
one dimensional lattice. In particular we proved that the 
transmission of a particle through an array of impuri- 
ties is independent of the direction of travel. This the- 
orem may help understand, among other things, weak 
localization, where time-reversed paths play an impor- 
tant role. Some important differences arise because of 
the lattice: first, it becomes clear that the group velocity 
is most important (in the continuum limit the energy is 
usually emphasized). Secondly, other symmetries exist 
on a lattice for particular wavevectors. Naturally these 
symmetries do not hold in the actual scattering region. 
While we have no specific proposals at present, we hope 
this work motivates experimentalists to look for these 
violations in the vicinity of particular impurity configu- 
rations. Further work is in progress with trimer impurity 
configurations, more general band structures, and higher 
dimensionality. 
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